Genetic differentiation and phylogeography of rotifer Polyarthra dolichoptera and P. vulgaris populations between Southeastern China and eastern North America: High intercontinental differences

Abstract Genetic differentiations and phylogeographical patterns of small organisms may be shaped by spatial isolation, environmental gradients, and gene flow. However, knowledge about genetic differentiation of rotifers at the intercontinental scale is still limited. Polyarthra dolichoptera and P. vulgaris are cosmopolitan rotifers that are tolerant to environmental changes, offering an excellent model to address the research gap. Here, we investigated the populations in Southeastern China and eastern North America and evaluated the phylogeographical patterns from their geographical range sizes, geographic–genetic distance relationships and their responses to spatial‐environmental factors. Using the mitochondrial cytochrome c oxidase subunit I gene as the DNA marker, we analyzed a total of 170 individuals. Our results showed that some putative cryptic species, also known as entities were widely distributed, but most of them were limited to single areas. The divergence of P. dolichoptera and P. vulgaris indicated that gene flow between continents was limited while that within each continent was stronger. Oceanographic barriers do affect the phylogeographic pattern of rotifers in continental waters and serve to maintain genetic diversity in nature. The genetic distance of P. dolichoptera and P. vulgaris populations showed significant positive correlation with geographic distance. This might be due to the combined effects of habitat heterogeneity, long‐distance colonization, and oceanographic barriers. Furthermore, at the intercontinental scale, spatial distance had a stronger influence than environmental variables on the genetic differentiations of both populations. Wind‐ and animal‐mediated transport and even historical events of continental plate tectonics are potential factors for phylogeography of cosmopolitan rotifers.


| INTRODUC TI ON
Two of the most important questions in ecology and biogeography are how and why species composition differs between geographic locations (Dambros et al., 2020). Environmental factors, geographic distance, and dispersal barriers are potential drivers. Some researches suggest that genetic diversity is related to geographical distances by a classical distance-decay relationship (Gómez et al., 2002;Gómez-Rodríguez et al., 2020). Others hold that niche differentiation including abiotic and biotic factors, not dispersal, controls the patterns of genetic differentiation (Jin et al., 2020;Li et al., 2020;Ma et al., 2019). But most studies have demonstrated that combining ecological and phylogeographic information is more effective in understanding the factors that drive the distribution of genetic diversity (Gabaldón et al., 2017). Cryptic species are difficult or even impossible to distinguish by their morphology. The putative cryptic species defined through phylogenetic models by clusters of sequences are called entities (Obertegger et al., 2015). Generalized mixed Yule coalescent (GMYC), automatic barcode gap discovery (ABGD), and Poisson tree processes (PTP) models are widely used approaches for cryptic entities delimitation (Fontaneto, 2014;Kordbacheh et al., 2017). Since the sequences are clustered by different models according to the genetic distance, they are called GMYC entities. (Obertegger et al., 2015). With the development of molecular tools, increasing numbers of cryptic entities have been discovered in various morphological species of rotifers (Fontaneto, 2014). Cryptic entities have been found in almost all groups of animals (Fossen et al., 2016;Tang et al., 2012), and rotifers seem to be one of the invertebrates hosting the highest potential cryptic diversity in the world (Fontaneto, 2014;Fontaneto et al., 2009). For instance, eight potential entities of Brachionus calyciflorus were found in eastern China (Xiang et al., 2011). Also, more than seven entities of Euchlanis dilatata were defined in North America (Kordbacheh et al., 2017). By the end of 2017, there were 15 entities of B. plicatilis recorded in the world (Mills et al., 2017). Rotifers are invertebrates sensitive to environmental changes. They are widely distributed all over the world and live in all kinds of water bodies (Liang et al., 2020). Thus, rotifers are advantaged for the study of phylogeography and excellent examples of cryptic entities. As the cosmopolitan and dominant species of rotifers show strong adaptability to the environmental changes, studying their spatial and temporal patterns of genetic differentiation is of great significance for understanding gene flow and adaptive evolution mechanisms in microscopic organisms (Zhang et al., 2018).
Rotifer biogeography is a complicated issue; biogeographical patterns exist but are difficult to detect and have been studied mostly at the continental scale (Fontaneto, 2019). To date, findings of the relationship between geographical distance and genetic distance are inconsistent. A study of Brachionus calyciflorus phylogeography in eastern China showed no significant association between geographical and genetic distances (Xiang et al., 2011). Moreover, a significant albeit weak correlation was found in Euchlanis dilatata phylogeography within North America (Kordbacheh et al., 2017). The relationship between geographical distance and genetic distance is far from clear in rotifers, as seen in studies of large geographic scales.
Baas-Becking's hypothesis, known as "everything is everywhere, but the environment selects" (EisE) posits that for small species, spatial variation occurs because of ecological differentiation and not because of restricted dispersal (Fenchel & Finlay, 2006). For example, the differential distribution of cryptic entities has been related to ecological factors: temperature and algal food concentration in the B. calyciflorus population (Li et al., 2010); salinity in the B. plicatilis population (Gómez et al., 2002); and total phosphorus concentration in the Synchaeta pectinata population .
Moreover, coexistence of cryptic entities is common in a single water body and even can be found in two distanlty separated locations Papakostas et al., 2016;Wen et al., 2016). However, most studies of environmental-spatial selection to date focused on small geographic or continental scales, and genetic differentiation across intercontinental scales has been underexplored. For instance, the two cryptic entities belonging to the Ascomorpha ovalis population occur at different altitudes in Mexico area (García-Morales & Elías-Gutiérrez, 2013). In another example, cryptic entities of a Polyarthra dolichoptera population was distributed along an altitudinal gradient in Trentino-South Tyrol, Italy.
Although several studies have supported the idea that most cryptic taxa within the population are widespread while others have locally restricted distributions (Fontaneto, Barraclough, et al., 2008;Mills et al., 2017), the knowledge about effects of geographical barriers on genetic differentiation of rotifer is still limited.
Polyarthra dolichoptera and P. vulgaris, cosmopolitan rotifer species, are more tolerant to seasonal changes than other rotifers and exist in a variety of freshwater bodies but not marine systems (Liang et al., 2019). They are zooplankton with limited mobility and can Wind-and animal-mediated transport and even historical events of continental plate tectonics are potential factors for phylogeography of cosmopolitan rotifers.

K E Y W O R D S
cryptic species, genetic distance, long-distance dispersal, phylogeography, Rotifer, spatiotemporal pattern

T A X O N O M Y C L A S S I F I C A T I O N
Population ecology be transported by currents within a single water body. As they are monogonont rotifers that reproduce by cyclical parthenogenesis, resting eggs (dormant stages) will be produced to escape unsuitable environmental conditions to temporarily avoid competition and parasitism and thus to move to distant areas and new habitats through abiotic and biotic factors (Fontaneto, 2019). It has been suggested that Polyarthra entities distribution reflects genotypic adaptations to temperature differences and food resources (Obertegger et al., 2015). Higher genetic diversity is seen in low-altitude lakes that are characterized by warm temperature and mesotrophic to eutrophic conditions (Obertegger et al., 2015). However, the major factors of genetic differentiation of rotifers at larger geographic scales is still unclear and needs further investigation.
In this study, we attempted to address whether oceanographic barriers affect the phylogeographic pattern and maintain biological diversity of rotifers in continental waters by investigating the genetic differentiation of P. dolichoptera and P. vulgaris populations from Southeastern China and eastern North America. The specific objectives were as follows: (1) to determine whether the populations of cosmopolitan species from different continents formed independent or overlapping strains; and (2) to understand whether environmental or spatial variables are the key factors affecting genetic differentiation.

| Sampling
Samples were collected from 28 sites including rivers, ponds, and lakes in both eastern North America and Southeastern China, during June 2018 to September 2019 ( Figure 1; Table 1). For determining detailed cryptic entities structure at a small geographic scale, five sites in Guangzhou City, China, and five ponds in New London Country, USA, were sampled. In consideration of temperature effects on entities structure, four seasons were sampled in Lake Liuye, China (6_liuye, 9_liuye, 12_liuye, and 3_liuye; Table 1).
All rotifer samples were collected by towing a plankton net (mesh size 30 μm) horizontally at surface and subsurface depths and preserved in 50-ml centrifuge tubes. Samples were fixed on site immediately with neutral Lugol's solution at 2% final concentration and transported in a cooler before storing at −20°C. In vivo semi-quantitative measurements of chlorophyll-a concentration (Chl-a) were obtained using a FluoroSenseTM handheld fluorometer (Turner Designs, USA). Water temperature (Temp) was measured on site. Total dissolved phosphorus (TP) and total dissolved nitrogen (TN) were determined in the laboratory following standard analytical methods (GB3838-2002, MEE, 2002. Also, GPS coordinates and altitude values were recorded.

| Species identification and isolation
Species identification was based on the most authoritative taxonomy review on Synchaetidae (Nogrady & Segers, 2002) using a regular compound microscope (Olympus BX51). Although P. dolichoptera and P. vulgaris are similar in morphology, they can still be distinguished by morphological features such as body sizes and fins. For P. dolichoptera, the main fins are as long or longer than body length; width/ length 1:10-1:15, serration coarse, widely spaced, uniformly swordlike with central rib almost reaching the end; serration lines faint or absent. For P. vulgaris, the main fins are as long or shorter than body length but broader than that of P. dolichoptera; width/length 1:5-1:7, feathery or leaf-like with strong central rib. Each identified individual was isolated with micropipette under the stereomicroscope (Olympus SZX16) ( Figure 2). Single individuals were rinsed several times and transferred into PCR tubes for DNA analysis.

| DNA extraction and amplification
DNA from each single animal was extracted following the HotSHOT protocol (Montero-Pau et al., 2008). Then, the partial mitochondrial cytochrome c oxidase subunit I (COI) mtDNA gene was amplified and sequenced using primers LCOI (5'-GGT CAA CAA ATC ATA AAG ATA TTGG-3') and HCOI (5'-TAA ACT TCA GGG TGA CCA AAA AAT CA-3') (Folmer et al., 1994). PCR was processed according to the TaKaRa exTaq protocol with 5 μl of extracted DNA. Cycle conditions were initial denaturation at 94°C for 3 min, followed by 35 cycles of denaturation at 94°C for 30 s, annealing at 52°C for 30 s, and extension at 72°C for 45 s. The amplification ended with a final extension of 72°C for 8 min. Successful amplification products were then purified using the TaKaRa Minibest agarose Gel DNA extraction Kit before being sent to TsingKe company for sequencing.

| Sequence alignment and phylogenetic analyses
Sequences were aligned in Mega X using Clustal-W and then checked manually (Kumar et al., 2018). The same program was used to translate COI gene sequences to proteins and to check for stop codons.

| Cryptic entities delimitation
An ultrametric tree generated by BEAST was required for both GMYC and PTP delimitations. The GMYC delimitation analysis was processed on software R 3.6.1 using the "rncl" and "splits" packages, (R Core Team, 2019). The GMYC model is a likelihood method for delimiting species by fitting within and between species branching models to reconstruct gene trees (Fujisawa & Barraclough, 2013;Pons et al., 2006).
The ABGD model was performed for primary species delimitation and was processed on the website https://bioin fo.mnhn.fr/ abi/publi c/abgd/ (Accessed January 10, 2020). ABGD classifies sequences into putative entities based on pairwise genetic distances without any prior assumptions (Puillandre et al., 2012).
PTP is a tree-based method that uses the number of substitutions to distinguish intraspecies processes from interspecies processes.
This method considers two classes of Poisson processes, speciation (higher substitution rate associated to interspecies events), and coalescence (within species events) (Zhang et al., 2013). PTP (http:// speci es.h-its.org/ptp/, Accessed January 8, 2020) was applied to ultrametric BEAST trees using the default settings to detect the number of entities (Kordbacheh et al., 2017).

| Geographical and genetic distance analysis
The geographic distance matrices were calculated in R 3.6.1 (R Core Team, 2019) using the "geosphere" package (Hijmans et al., 2019). In order to explain the relationships between geographic distances and genetic distances, linear regressions were processed in R 3.6.1. version 1.2.3 of the Geographic Distance Matrix Generator was used to construct a geographic distance matrix. Mantel tests were run in the R package "ecodist" with 10,000 permutations to test whether genetic variation among populations is correlated to geographic distances among populations (Kordbacheh et al., 2017). To determine the significance of differences (p < .05) in genetic distances among three different F I G U R E 2 Morphologic photographs: P. vulgaris (a); the trophi of P. vulgaris (b); P. dolichoptera (c); the trophi of P. dolichoptera (d) groups, analysis of variance (ANOVA) with TukeyHSD test was conducted, using the software R 3.6.1, "agricolae," "car," and "multcomp" packages.

| Relationships between entities and environmental factors
Since the approximation by the linear function is poor over a long gradient when the data are too heterogeneous, detrended correspondence analysis (DCA) should be carried out. CCA or RDA model is determined based on the 27 sample datasets (Table 1)

| Genetic diversity
We obtained 107 COI sequences of P. dolichoptera and 63 sequences of P. vulgaris, with aligned lengths of 562p and 589 bp, respectively (Accession numbers:

| Phylogenetic analysis and entities delimitation
Using Bayesian phylogenetic analysis, P. dolichoptera population was reconstructed ( Figure 3) A large number of cryptic entities were detected using three independent methods (Figures 3 and 4). The ABGD method produced 20 cryptic entities for P. dolichoptera population and four for P. vulgaris. GMYC analysis revealed 21 cryptic entities in P. dolichoptera population and five in P. vulgaris population. Using the PTP method, P. dolichoptera population was delimited into 24 entities and P. vulgaris population into 13. The most conservative estimate of cryptic entities was obtained using ABGD, while the PTP method gave the greatest number of cryptic entities. All three methods shared common species boundaries for the smallest number of cryptic entities, however. As the results from ABGD and GMYC were similar, unless specified otherwise, cryptic entities determination will be discussed based on ABGD results (ABGD entities).

| Geographical distribution and coexistence
The range size of geographic distances declined as the resolution of classification increased from morphological species to cryptic entities to haplotype ( Figure S1). Both P. dolichoptera and P. vulgaris species were widely distributed, with the geographic range sizes up to 12,903 and 12,836 km, respectively. Their range sizes at cryptic entities and haplotype levels decreased to 2726 and 411 km, respectively. These correspond to the mean range sizes at species level (P. dolichoptera: 5374 ± 5665 km; P. vulgaris: 5517 ± 5720 km) significantly larger than at entities (463 ± 815 km) and haplotype level (12 ± 61 km) (p < .05). However, there was no Coexistence and seasonal succession also occurred in the present study. D1 and D15 coexisted in June in Liuye Lake, while D1, D3, and D4 coexisted in September in the Chuanzi River (Figure 3). V1 and V2 also coexisted in June in Liuye Lake (Figure 4). On the other hand, D8 only appeared in December in the Chuanzihe River.

| Phylogeographical patterns and genetic structure
The relationships between dependent variables for genetic distance and independent variables for geographic distance were examined by linear regression. The genetic distance of P. dolichoptera population showed significant positive correlation with geographic distance (R 2 = .18, p < .01) (Figure 5a). The genetic distance of P. vulgaris population was also positively correlative with geographic distance (R 2 = .53, p < .01) (Figure 5b).
Our results indicated that in both P. dolichoptera and P. vulgaris populations, the mean genetic distances between the two continents were significantly higher than those within continents were significantly higher than those within a single continent (p < .05). This suggests that there is a higher level of gene flow and higher frequency of recombination within continents than between continents.

| Relationships among entities distributions, spatial and environmental factors
As the longest gradient performed by detrended correspondence analysis (DCA) was 7.5 (larger than 4, Table S2), a canonical correlation analysis (CCA) model was chosen for estimating the relationship among entities (within the P. dolichoptera and P. vulgaris), spatial and environmental factors. The first two ordinate axes explained 49.5% of the entities distribution, spatial and environment variability in the CCA ordination (Table S3). The correlation coefficient of CCA ordination showed that three variables including longitude, latitude, and altitude were significantly related to the entities distribution (p < .05) (Table S4). However, the environmental factors that include temperature, chlorophyll-a, TN, and TP concentration were not significant variables affecting the entities' distributions (Table   S4). Figure 7 clearly showed that the Southeastern China samples (red) mostly stayed on the left of the figure, while the eastern North America samples (blue) were mostly on the right of the figure. In addition, the spatial variables of longitude were the longest and showed positive correlation with axis 1 (with the 95% eigenvalues), which indicated that longitude was the key factor for the variation of the entities' distributions for both species. Furthermore, the entities-variables relationship was similar to that of sampling sitesvariables, which indicated that most entities (within the P. dolichoptera and P. vulgaris) tended to be restricted to specific regions.

| Selection of DNA markers for genetic differentiation studies
Nuclear (18S, 28S) ribosomal RNA genes, nuclear internal transcribed spacer (ITS), and mitochondrial cytochrome c oxidase subunit I gene (COI) are widely used as DNA barcodes for identification (Papakostas et al., 2016). COI has been used as the most popular marker of genetic diversity in animals and was elected as F I G U R E 5 Linear regression model between the uncorrected pairwise genetic distances and geographic distance of the P. dolichoptera populations (a) (n = 5724) and the P. vulgaris populations (b) (n = 1984) the standardized tool for molecular taxonomy and identification over the last three decades (Ratnasingham & Hebert, 2007). In addition, COI showed great applicability in the study of zooplankton gene flow and thermal adaptation (Sasaki & Dam, 2019). Our results also showed that COI performed well in genetic diversity studies of P. dolichoptera and P. vulgaris.
Using one single marker, especially COI usually causes some problems due to its vulnerability to the selection process. The lack of recombination for this marker makes the data sensitive to selective sweeps. This is not ideal if one wishes to investigate the demographic processes in a population (Galtier et al., 2009). Also, the use of 18S rDNA for species in rotifers is discouraged, because of its low variability (Tang et al., 2012). ITS has been suggested to be a more reliable marker for cryptic species delimitation of Brachionus (Mills et al., 2017;Papakostas et al., 2016). But given that COI is more variable than ITS, the former is still the best marker to be used for exploration of population genetic structure within species and phylogeography (Mills et al., 2017). So it is better to combine COI with internal transcribed spacer (ITS) for phylogeography and cryptic species delimitation (Kordbacheh et al., 2017).

| Genetic divergence and geographic distribution
Long-distance dispersal of entities has been reported not only in Monogononta including Brachionus, Polyarthra, Euchlanis, and Lecane but also in Bdelloidea including Philodina and Rotaria Kordbacheh et al., 2017). The entity D1 was found in Guangdong and Hunan provinces, separated by >500 km (e.g., liuye,

F I G U R E 6
The ranges of uncorrected pairwise genetic distances among the three groups of P. dolichoptera (a) and P. vulgaris (b). Differences were detected with the Turkey HSD method. Letters indicate sample means that are similar (same letter) or significantly different (different letter)

F I G U R E 7
Canonical correlation analysis (CCA) of the P. dolichoptera and P. vulgaris entities with environmental and spatial variables. Relationship between sampling sites and variables (a). Relationship between entities and variables (b). Red, samples and entities from Southeastern China; Blue, samples and entities from eastern North America. Abbreviations used in the figures: V, P. vulgaris entities; D, P. dolichoptera entities; 3_liuye, in Lake Liuye in March; 6_liuye, in Lake Liuye in June; 9_liuye, in Lake Liuye in September; 12_liuye, in Lake Liuye in December; 9_chuanzi2, in the Chuanzi River in September; 12_chuanzi2, in the Chuanzi River in December chuanzi2, minghu, pearlriver2). Also, D10 was found across a range of >500 km in the US states of New York and Connecticut (e.g., SUC, Niagara, P2, P4). In addition, V1 was widely distributed in southeastern of China, while V3 was widespread in Connecticut.
The cosmopolitan distribution of rotifers could be attributed to long-distance dispersal. Colonization and long-distance dispersal to different waters across whole or even multiple continents, which may be mediated by waterfowl, have been observed in a number of zooplankton species (Gómez et al., 2002). Some areas that are widely separated share haplotypes and therefore appear genetically connected (Xiang et al., 2011). Sasaki and Dam (2019) found that the widely distributed genetic clades of a marine copepod shared haplotypes between geographically distant populations. This implies that gene flow can be strong enough to overcome long distances at least within a continent.
Most small organisms do have very widespread distributions, but some are limited to distinct geographical areas (Fontaneto, Barraclough, et al., 2008;Savary et al., 2018 (Xiang et al., 2011). Although high genetic distances of Adineta can be found at different geographical distances, closely related individuals were only found at geographical scales <2000 km (Fontaneto, Barraclough, et al., 2008). In the current study, we found that the range sizes of geographic distance in both genera declined as resolution increased from species to entities to haplotype. This suggests that the restricted entities in our study are not simply an artifact of sampling fewer individuals at lower levels.
Interestingly, even though Polyarthra was widely distributed as a genus, no identical entities or haplotypes have been found to appear on both continents at the same time. High levels of nucleotide diversity could be caused by a particularly high mutation rate which results in genetic divergence (Xiang et al., 2011). The nucleotide diversity of the total P. dolichoptera and total P. vulgaris samples were relatively higher than those from the single continent, which suggested high genetic divergence between the samples from the different continents. Given that the small sample sizes in this study impeded a thorough detection of genetic diversity, the degree of cryptic diversity in Southeastern China and eastern North America is likely to be higher than what we reported here. This genetic diversity could explain why cosmopolitan rotifers have apparently wide tolerance toward spatial and environmental changes.
Our results indicated that all of the entities from eastern North America formed independent strains that were separate from the Chinese ones, indicating high divergence. These results are consistent with the study in Noctiluca, a heterotrophic dinoflagellate (Pan et al., 2016). The haplotypes of Noctiluca within China were geographically quite homogeneous, but were generally different, compared to the American population, suggesting basin or continental-scale endemism. In addition, a study of the entities of B. plicatilis revealed existence of four clades associated to four geographic regions (one in North America, two in Europe and one in Australia) (Mills et al., 2017).
Levels of gene flow can be estimated by producing visible patterns using allele frequencies and DNA sequence differences (Slatkin, 1987). Lack of differentiation in mitochondrial COI sequences of geographically distant populations usually indicates strong effects of gene flow (Sasaki & Dam, 2019). Our study showed that the genetic distances of Southeastern China VS eastern North America were significantly higher than those within each continent.

| Relationship between geographic and genetic distance
In the present study, a significantly positive correlation between genetic and geographic distance was found in both Polyarthra species.
Similar results have been obtained for E. dilatata in North America (Kordbacheh et al., 2017). Moreover, it was reported that there was a strong positive correlation between genetic distance and geographic range when comparing samples at small geographical scales (Kordbacheh et al., 2017). Habitat heterogeneity and temporal variation can generate high genetic diversity on small geographic scale study (Fontaneto et al., 2009).
However, as the geographical scope of the study becomes broader, this correlation may weaken or disappear. In another example, no significant associations between geographical and genetic distances were found for B. calyciflorus across eastern China.
The nonsignificant correlation may result from the effects of longdistance colonization and secondary contact, combined with monopolization effects which reduce gene flow among established populations (Kordbacheh et al., 2017;Xiang et al., 2011). In our results, Guangdong and Hunan provinces shared P. dolichoptera entities D1 and P. vulgaris entities V1, while New York and Connecticut states shared the D10 and V3. Thus, long-distance intra-continental dispersal and colonization are responsible for depressing the geographic-genetic correlation in the present study.
The studies of phylogeography of rotifers usually show a weak correlation between geographical and genetic distances within continents (Kordbacheh et al., 2017;Xiang et al., 2011). Since the between-continent genetic distances were significantly higher than those within continent, a stronger positive correlation between genetic and geographic distance was observed, in consideration of datasets from the two continents. Furthermore, the significant genetic difference between the trans-Pacific regions suggests gene flow limitation (Pan et al., 2016). Therefore, effects of the Pacific barrier leads to the restriction of gene flow and results in an increase in genetic distance.
The relationship between geographic distance and genetic distance may not be simply linear. The results of linear regression showed that genetic distance increases rapidly with geographical distance at small geographic scales because of the coexistence and habitat heterogeneity among lakes, streams, and rivers. As the geographic distance extends to an entire continent, long-distance colonization leads to the decrease of the genetic distance. Thus, a power-law model is expected when there is no dispersal limitation (Gómez-Rodríguez et al., 2020). However, extreme barriers to dispersal, such as separate oceanographic basins, lead to an increase in genetic distance at larger geographic scales. Outliers above the correlation line between geographic distance and genetic distance might suggest a significant dispersal barrier at large geographic scales.
However, the small sample sizes in this study impeded a definitive evaluation of genetic diversity. A more complete picture of the genetic and phylogeographic patterns of the rotifers in question could be established if samples were also obtained on Africa, Oceania, and the west coast of North America.

| Key factors for phylogeographical patterns of rotifers
CCA indicated that spatial variables including longitude, latitude, and altitude were key factors in controlling the entities' structure rather than environmental factors such as temperature TP, TN, and chlorophyll-a concentrations. For most populations, genetic differentiation depend largely upon the evolutionary force regulating spatial patterns rather than seasonal differentiation (Obertegger et al., 2015;Xiang et al., 2011). In one study, cryptic diversity of P.
dolichoptera was found to change along an altitudinal gradient in the Trentino-South Tyrol region, but environmental parameters such as temperature and trophic status might also affect the distribution of entities (Obertegger et al., 2015). Thus, in the absence of geographical barriers, genetic divergence might be more explained by environmental gradients (Tisthammer et al., 2020).
The present study compared populations with similar food resource levels (Chl-a) at different sites and we found that the haplotypes belonged to different clades. Therefore, food sources level might not be an influencing factor for their genetic divergence.
Dispersal, genetic diversity and gene flow can be strongly affected by temperature changes (Sasaki & Dam, 2019). In our study, no winter samples were collected because lakes in eastern North America were frozen over the winter. Although the ambient temperatures of our samples ranged from 8 to 30℃ (Table 1), Polyarthra rotifers in eastern North America actually experienced even greater interannual changes in temperature. In the future, high-frequency and long-period sampling surveys can be carried out in a single water body to further explore the influencing factors of the coexistence and seasonal succession of Polyarthra spp.
Although freshwater rotifers are already confined to an enclosed area, gene flow can be increased in geographically distant areas within the continent through wind-, water-, and animal-mediated factors. Nevertheless, a quantification of the role of such potential long-distance dispersal on the biogeography of microscopic aquatic animals is unavailable now (Fontaneto, 2019). Rotifers possess the ability for passive long-distance dispersal through their diapausing stages including diapausing eggs and xerosomes .
Even a small number of successful dormant propagules is predicted to allow for effective long-distance dispersal, especially if those propagules can rapidly develop large populations via parthenogenesis (Fontaneto, 2019). In this way, hydrology influencing community composition and wind influencing dispersal could also play an important role in rotifer dispersal (Liang et al., 2019;Rivas et al., 2018).
In addition, as boating is one of the popular recreational activities for Americans in summer, rotifers and resting eggs can spread over North American lakes by launching boats. Human-mediated transport has likely facilitated the species' persistence since its initial colonization, through the ongoing introduction and intercontinental spread of genetic variation (Baird et al., 2020). Wind-and migratorymediated transport could operate at larger scales but are impeded by the oceanic barriers. Since dust storms and waterfowl cannot cross the strong oceanographic barriers, gene flow via dispersal of diapausing eggs and xerosomes among rotifer populations could be limited.
The phylogeographic pattern of Polyarhtra spp. at a global scale may also be related to the historical events of continental plate tectonics. The biogeographic realms of rotifers can be defined into eight regions (e.g., Palaearctic, Nearctic, and Oriental) (Segers, 2007). During the Cretaceous, the North American continent was still connected to the Eurasian continent, and zooplankton such as rotifers spread widely in Laurasia (Segers, 2007;Wu, 2012). In the early Cenozoic, as the North American continent drifted away from the Eurasian continent, gene flow in the established populations between the two continents was likely impeded. Southeastern China, now part of the Oriental Region, originally belonged to Gondwana. It was closely integrated with the Eurasian continent after the collision of the Indian subcontinent. Subsequently, a large amount of long-distance dispersal and gene flow took place between the rotifer species in the Oriental region and the Palaearctic region (Wu, 2012). Since the time scales involved in the movement of the continents via plate tectonics are extremely long, however, contemporary distribution are more likely due to recent dispersal mechanisms.

| CON CLUS IONS
1. Cryptic diversities of Polyarthra dolichoptera and P. vulgaris are definitely underestimated in the world. Oceanographic barriers do affect the phylogeographic pattern of rotifers in continental waters and serve to maintain genetic diversity in nature.

Entities of P. dolichoptera and P. vulgaris from eastern North
America formed independent strains that were separate from the Chinese ones, indicating high divergence. The divergence indicates that gene flow between eastern North America and Southeastern China is limited while that within eastern North America or Southeastern China was higher.
3. The genetic distance of both P. dolichoptera and P. vulgaris populations showed significant positive correlation with geographic distance at the intercontinental scale. This may result from the effects of habitat heterogeneity, long-distance colonization, and oceanographic barriers to dispersal.
4. Spatial variables are key factors in affecting the genetic differentiation of rotifers when compared with environmental variables at the intercontinental scale. Wind-and animal-mediated transport and even historical events of continental plate tectonic are potential factors for phylogeographic patterns of cosmopolitan rotifers.

ACK N OWLED G M ENTS
We thank the National Natural Science Foundation of China Xiaotong Ye for sampling assistance and sequencing.

CO N FLI C T O F I NTE R E S T
None declared. Yufeng Yang: Project administration (lead); Supervision (lead).

DATA AVA I L A B I L I T Y S TAT E M E N T
The COI sequences lacked internal stop codons of this study have been submitted to the NCBI database. The GenBank accession numbers are listed in Table S1.